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Abstract. We show that the dynamics of the surface plasmon in metallic nanoparticles damped by its 
interaction with particle-hole excitations can be modelled by a single degree of freedom coupled to an 
environment. In this approach, the fast decrease of the dipole matrix elements that couple the plasmon to 
particle-hole pairs with the energy of the excitation allows a separation of the Hilbert space into low- and 
high-energy subspaces at a characteristic energy that we estimate. A picture of the spectrum consisting 
of a collective excitation built from low-energy excitations which interacts with high-energy particle-hole 
states can be formalised. The high-energy excitations yield an approximate description of a dissipative 
environment (or "bath") within a finite confined system. Estimates for the relevant timescales establish the 
Markovian character of the bath dynamics with respect to the surface plasmon evolution for nanoparticles 
with a radius larger than about 1 nm. 

PACS. 73.20.Mf Collective excitations - 78.67.n Optical properties of low-dimensional, mesoscopic, and 
nanoscale materials and structures - 71.45.Gm Exchange, correlation, dielectric and magnetic response 
functions, plasmons 



1 Introduction 

One of the main questions driving the experimental and 
theoretical research of the last twenty years on metallic 
nanoparticles is how large in size do we have to go to ob- 
serve bulk properties. The answer of course depends on 
the physical property under study. It is furthermore never 
clearcut. But in many cases the size needed to observe 
bulk- like behaviour is surprisingly small. Furthermore, the 
size dependence of quantitative features can already be 
smooth for very small systems. For instance, the binding 
energy of Nag clusters is quite close to that of Nas , even 
if these two systems are very different from the molecu- 
lar point of view [1]. Such a continuity points towards the 
relatively minor importance of the ionic cores and sup- 
ports the descriptions based on the jellium model, where 
the conduction electrons are subject to a uniform neu- 
tralising background [2] . This view has been validated by 
the evidence of electronic shells provided by the abun- 
dance spectra of alkaline clusters [3]. Within the jellium 
approximation, the density of electronic states is given by 
a bulk-like contribution to which we have to add surface 
and periodic-orbit (or shell) corrections [4]. 

The optical properties of metallic clusters are domi- 
nated by the response of the conduction electrons. Ex- 



cept for the smallest clusters, where the transitions be- 
tween single particle levels dictate the optical response, 
the optical absorption is dominated by a collective ex- 
citation, the surface plasmon. This resonance is located 
near the classical Mie frequency ojm — w p /VI, where 
uj p = (47rn c e 2 /TO c ) 1 / 2 is the bulk plasma frequency and e, 
m e , and n c denote the electronic charge, mass, and den- 
sity, respectively. Nevertheless, surface effects lead to a 
reduction of the surface plasmon frequency with respect 
to the bulk value ujm [2,5,6]. The question of how large 
the size of the nanoparticle has to be in order to observe 
a collective excitation in a finite system has the surpris- 
ing answer that a very small cluster, like for instance Na@, 
may already be enough [7,8]. 

For nanoparticle radii a between about 0.5 and 5nm 
(N = 8 to 14000 conduction electrons for the case of Na) 
the main effect limiting the lifetime of the resonance is the 
Landau damping, i.e., decay into particle-hole pairs. The 
resulting linewidth is given by 



7tot(o) -7i + 7(«)+7 osc (a) 



(1) 



where ji is the intrinsic bulk-like linewidth. The second 
term in the right-hand side of fl} decreases with the na- 
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noparticle size as [9-11] 



3u F 

7(a) = 1^° 



M 



(2) 



with £p = h 2 kp/2m c and up the Fermi energy and veloc- 
ity, respectively. The monotonically increasing function go 
is defined as 

9o(x) = j^{^x{x + l)(4x(x + 1) + 3) 
- 3(2x + 1) In (v 7 ^ + Vx + T) 



\/x(x- l)(4x(x - 1) + 3) 
3(2x- 1) In (a/x + v / ^ 3 T) 6»(x- 1)} , (3) 



being the Heaviside step function. The last contri- 
bution 7 osc ~ cos(2fc F a)/(fc F a) 5/2 to the linewidth (JU) is 
nonmonotonic in a, and arises from the correlation of the 
densities of states of the particles and holes [12,13]. It 
becomes relevant for the smallest sizes of the considered 
interval, typically for radii in the range 0.5-1.5nm. For 
radius larger than 5nm the Landau damping competes 
with radiation damping in limiting the lifetime of the res- 
onance. For radii smaller than 0.5 nm the interaction with 
the ionic background becomes dominant, and the jellium 
model no longer provides a useful description. 

The use of femtosecond pulsed lasers in pump-probe 
spectroscopy has rendered possible to experimentally ad- 
dress the surface plasmon dynamics [14-17]. The theoret- 
ical descriptions built to study this problem [5, 6, 13, 18] 
treat the collective coordinate as a special degree of free- 
dom, which is coupled to an environment constituted by 
the degrees of freedom of the relative coordinates. Fric- 
tion arising from particle-hole pairs has been studied in 
bulk systems [19]. In the present work we are interested 
in the case where the excitation spectrum arises from a 
finite number of particles. The finite number of electrons 
makes the term "environment" not completely justified 
in this situation, and leads us to reformulate the above- 
mentioned question in dynamical terms as: How large in 
size do we need to go to be allowed to describe the rel- 
ative coordinates of the electron gas as an environment 
damping the collective excitation? 

This last question is the main subject of the present 
paper. In particular, we determine which are the energies 
of the electronic excitations that are active in the damp- 
ing of the surface plasmon, and which is the characteris- 
tic response time of the large enough electronic environ- 
ment. The latter is important in justifying the Markovian 
approximation that is assumed to describe the dynam- 
ics of the surface plasmon coupled to particle-hole excita- 
tions [20]. 

The paper is organised as follows: In Section [5] we re- 
call the random phase approximation for the surface plas- 
mon and the separation of the excitation spectrum in low- 
and high-energy particle-hole excitations. In Scction[3l we 
show how the surface plasmon is built from low-energy 
particle-hole excitations and in Section H] the cutoff energy 



separating the above-mentioned low- and high-energy sec- 
tors is determined. In Section [5] the time dynamics of the 
environment giving rise to the damping of the collective 
excitation is considered. We end with the conclusions in 
Section 



2 Random phase approximation for the 
surface plasmon 

The understanding of collective excitations in metallic clus- 
ters has benefited from the accumulated knowledge in the 
related problem of nuclear physics, i.e., the description of 
the giant dipolar resonance [21,22]. The simplest many- 
body approach yielding collective excitations in a finite 
system is the Tamm-Dancoff description, where the ex- 
cited states are built from the Hartree-Fock ground state 
and all possible one particle-one hole (lp-lh) excitations. 
The basis considered is then mixed by the matrix elements 
of the residual interaction 



t t 



Vhf, 



(4) 



where cj, (c a ) creates (annihilates) the single-particle state 
| a) of the Hartree-Fock problem. In the above equation, 
v a 0-f5 = Va(3-y8 — Va(3S-y with w Q /3 7 5 the two-body matrix 
element of the Coulomb interaction, while 



1 



VW = - 



4 ^ 

a/3 



(5) 



is the Hartree-Fock interaction Hamiltonian. Diagonaliza- 
tion within this reduced Hilbert space gives the excitation 
energies. This task is considerably simplified under the as- 
sumption that the matrix elements v a p-yS of the residual 
interaction ^ are separable [23], 

(6) 



where d a ~ t = (ev|z|7) are dipole matrix elements and A is a 
positive constant characterising the repulsive residual in- 
teraction. Within the Tamm-Dancoff approximation, the 
excitation energies are given by the secular equation 



1 
A 



E 



E- As 



ph 



(7) 



where p (h) denotes a particle (hole) state above (below) 
the Fermi level, and Ae p h = e v — £h are particle-hole (p-h) 
excitation energies. 

One important drawback of the Tamm-Dancoff ap- 
proximation is that its ground state does not include cor- 
relations. This is not the case with the random phase ap- 
proximation (RPA) where the basis of the reduced Hilbert 
space is built from all possible lp-lh creations and de- 
structions acting on the ground state. Using again the 
separability hypothesis ^ for the residual interaction, the 
RPA secular equation 



1 

A 



S(E) 



(8) 
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is different from ([7]). Here we have denned the RPA sum 

2Ae p h\d p h\ 2 



Ac 

in 1 1 l /e F 



4-10 



E 2 - Ae 2 h 
P h P h 



(9) 



The validity of the RPA [8] and of the separability hypoth- 
esis within the time-dependent local density approxima- 
tion [24] have been well established, even in small clusters. 

Nuclear physics textbooks (see, e.g., Ref. [21]) show 
how to solve the secular equations ([7]) and graphi- 
cally under the implicit assumption that the p-h spec- 
trum is bounded in energy. In that case, the solutions of 
the secular equations are merely renormalizations of the 
p-h energies Ae p h except for the largest excitation energy 
which corresponds to the collective mode. However, it is 
clear that in a metallic cluster the p-h excitations are not 
bounded from above on the scale of the plasma energy. As 
we will see in the sequel, it is the fast decay of the dipole 
matrix elements d p t with the p-h energy that ensures the 
applicability of the standard picture. 

An ingenious way of describing the collective excita- 
tion in metallic clusters which circumvents the problem 
of the unbound p-h spectrum is the separation of the re- 
duced lp-lh RPA Hilbert space in a low-energy sector 
(the restricted subspace), containing p-h excitations with 
low energy, and a high-energy sector (the additional sub- 
space) [11,22]. The residual interaction gives rise to the 
collective excitation as a coherent superposition of a large 
number of basis states of the restricted subspace. This 
excitation energy lies in the high-energy sector, and the 
nonvanishing coupling with p-h states of the additional 
RPA subspace results in the broadening of the collective 
resonance. In this approach the cutoff energy separating 
the two subspaces is somehow arbitrary, as long as it is 
taken smaller than the resonance energy. The arbitrari- 
ness of the cutoff is not problematic for calculating the 
lifetime of the resonance [11], but the approximate theo- 
retical description of other physical quantities might de- 
pend on the cutoff. For example, in the approach where 
the environment is separated from the low-energy excita- 
tions that compose the collective plasmon excitation, the 
timescales that characterise the dynamics of the electronic 
environment [20] depend on the cutoff energy. Moreover, 
it was shown in [6] that the environment-induced redshift 
of the resonance frequency depends logarithmically on the 
cutoff. In the following we will provide an estimation of 
the cutoff energy separating the two subspaces. 




Fig. 1. RPA sum (|13|1 for a Na nanoparticle with k^a = 30 
(solid black line). The p-h excitation energies Ae p h (repre- 
sented by vertical grey lines) have been obtained from the 
semiclassical spectrum (|4ip . There are about 15000 degener- 
ate excitations in the interval shown, and the smallest excita- 
tion zAe m m ~ £f/10 has a degeneracy M = 380. The dashed 
line takes only into account the contribution of Ziemin (see 
Eq. fl24j l). The cutoff energy Ae c separating the two RPA sub- 
spaces is also shown in the figure (see Eq. (|30[l ). Above Ae c the 
two sums are close to each other, showing that for high energies 
the RPA sum is essentially given by the contributions coming 
from the low-energy p-h excitations. The horizontal dotted line 
indicates the position of the coupling constant A entering the 
RPA secular equation (J5J), according to the estimate (|25|l . 



The angular part is expressed in terms of Wigner-3j sym- 
bols as 



A 



m p 7nh 
l l 

L p L h 



(-l) m *J(2l p + l)(2l h + l) 



ft 


h 


i 




1 









\-m p m h 





(11) 



and sets the selection rules l p = 1% ± l and m p — rrih for 
the total and azimuthal angular momenta, respectively. 
The radial part depends on the energies of the p-h states 
as 

ni ^ Sh) = Wa^- (12) 

With the help of equations fTU)) - (fT2")) and using the 
appropriate dipole selection rules, we can write the RPA 
sum ((9|) as 



3 The surface plasmon as a superposition of 
low-energy particle-hole excitations 

The secular equations ([7J and ([5]) depend crucially on the 
form of the dipole matrix element d p h ■ Assuming that the 
p-h states are confined within a hard-wall sphere, one can 
decompose [11] 



Lph 



m p m h 



p 1 1 1 



(10) 



S(E) = 



,l p — Zfcil 



{E 2 - As 2 



ph 



where and n p are radial quantum numbers. We have 
defined fi p = h + 1 if l p = h + 1 and ft = l h if l p = l h - 1. 
Using the semiclassical quantisations (14TH and (|4ip of Ap- 
pendix El we finally obtain the result shown in Figure [1] 
In the figure we show (solid black line) the resulting 
i?-dependence of (fT5|) . The crossings of this curve with 
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the horizontal dotted one with height 1/A yield the exci- 
tation spectrum within the RPA. The lowest p-h energy 
estimated in Appendix [X] is 



Ae r 



£f 



kpa/ir 



(14) 



Whenever the energy E coincides with a p-h excitation en- 
ergy Ae v h, we have a divergence in S{E) (see vertical grey 
lines in Fig. [I]). For the lowest energies E of the interval 
considered the sum is dominated by the term associated 
with the divergence closest to E. On the other hand, for 
the largest energies E, the fast decay of d p h with Ae p h 
means that the divergences are only relevant for energies 
extremely close to them. Away from the divergences the 
sum is dominated by the contributions arising from the 
low-energy p-h excitations. 

To gain physical insight and proceed further in the 
analysis of the electron dynamics in metallic nanoparti- 
cles, we introduce the typical dipole matrix clement for 
states separated by a given energy difference As, 



d?- h (Ae) 



ph 



1/2 



pP- h (Ae) J ep _ A 



-1 1/2 



deC(e,Ae) 



(15) 



Here, we have introduced the local density of dipole matrix 
elements 

C(e, Ae) = \d P h\ 2 6(e ~ e h )6(e + Ae - e p ) . (16) 

ph 

As shown in Appendix [Bj this can be expressed as 



^ = 3^(4). 



(17) 



where 



F(x) = (2x + l)y/x(x + l) - ln(V^+ VxTT) . (18) 
For Ae -C e, equation (fT7]) simplifies to 

2a 2 e 2 



C(e,Ae) 



3tt 2 Ae 4 ' 



(19) 



In expression (fl"5 



pP- h (Zk) = ]T S(Ae ph - Ae)S lh , lp±1 S mh , mp (20) 

ph 

is the density of p-h excitations with energy Ae respect- 
ing the dipole selection rules. In order to simplify the pre- 
sentation, we do not consider spin degeneracy factors. In 
Appendix [Cl we show that for Ae <C ep we have 



p^{Ae) 



(hp a) Ae 



(21) 



and therefore in such a limit the typical matrix element 
()15l) can be approximated by 



dP- h {Ae 



^3k F a \Ae) 



(22) 



As a check for our estimation of the typical dipole matrix 
element, we have evaluated the energy- weighted sum rule 
J2 P h Ae p h\d p h\ 2 using (|2"2"1) and obtained about 70% of the 
exact result (3/4ir)ti 2 N/2m e [2]. This is quite reasonable 
regarding all the approximations we made to obtain (|22p . 

In order to illustrate the importance of the low-energy 
p-h excitations, we present in Figure [1] (dashed line) the 
contribution to the RPA sum coming only from the in- 
frared p-h excitation energy with the appropriate degen- 
eracy factor J\f (see Eq. (|43|) in App. IA"|) . Indeed, we can 
estimate © as 



S{E) m Af 



2Ae min [dP- h (Ae min )Y 
E 2 - AeL„ 



With the results flijl. ((2"2"j) . and we obtain 



S(E) 



E 2 - (e F 7r/fc F a) 2 



(23) 



(24) 



While in the lower part of the energy interval in Figure [T] 
the two curves exhibit considerable discrepancies, in the 
second half of the interval they are very close (except of 
course at the divergences). Since the collective excitation 
is found in this last interval, we see that it is mainly the 
low-energy p-h excitations that are relevant for the defi- 
nition of the collective excitation. 

Since the resonance energy is known from experiments, 
we can obtain the value of the coupling constant A. Indeed, 
using the estimate ([2~4l evaluated at E = Itwm in the sec- 
ular equation ffSj) , we obtain 



UK 



3tt 5 



-(k-pa) 



(Sw m ) 2 



(25) 



to leading order in kpa ^> 1. Notice that this result is 
consistent with the estimate obtained from the energy- 
weighted sum rule given in [11]. For the case studied in 
Figure [U our estimate ([2"B"j) yields k F ep/\ » 1600. As the 
radius a increases the lowest p-h energy decreases like 1/a, 
the degeneracy at this value increases as a 2 , and the den- 
sity of p-h excitations contributing to (JSJ, p p ~ h (Ae), grows 
as a 4 . This increase in the number of excitations contribut- 
ing to the sum is partially cancelled by the 1/a behaviour 
of the typical dipole matrix element, resulting in a de- 
crease of the coupling constant A proportional to the num- 
ber of particles in the cluster, A ~ 1/a 3 , and a value of 
the plasmon frequency which remains almost unaffected. 
Therefore, for a larger nanoparticle size, the divergences 
shown in Figure [T] would be more dense in energy, starting 
at a lower energy, and the vertical scale would increase as 
a 3 , thus obtaining a similar value for the plasmon excita- 
tion energy Hum- 
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4 Separation of the reduced and additional 
particle-hole subspaces 

As we have shown in the last section, the high-energy part 
of the p-h spectrum is not crucial for the determination 
of the energy of the collective excitation. In what follows 
we make this statement more quantitative and estimate 
the upper-bound cutoff Ae c of the low-energy excitations 
that we need in order to obtain a stable position of the 
surface plasmon. 

In order to obtain a quantitative estimate of the cutoff, 
we require that by changing it from Ae c to (3/2)Z\e c , the 
position of the plasmon changes only by a fraction of its 
lincwidth 7, the smallest energy scale with experimental 
significance. Our criterion leads to the condition 



Sac c (Swm) = Sa A (fiw M + fry) 



(26) 



with the RPA sum S that has been denned in (|9|). The 
additional subscript refers to the upper bound of the p-h 
energies. 

The left-hand side of (|2"fH) can be estimated according 

to 



dAep^(Ae) [ 1 >[ 
(Awm) — Ae z 



(27) 

with d p ~ h and p p ~ h as defined in (115j) and (|2"U|) , respectively. 
Using (T3TJ> and ([22]) . we obtain 



4fc 



' k F a 



1 



1 



37T 2 V hjjJM ) V ^ £ min As, 

to leading order in fcpo. Similarly, 



(28) 



in 1 \ TujmJ V w m 
1 2 

x 



As n 1 1 1 



3Ae ( 



(29) 



Using the expressions © and (fl4"|) for 7 and Ae m i n , re- 
spectively, finally yields according to the criterion ([26]) the 
cutoff energy 



Ae c 



9g (e F /ftujM) 



M : 



(30) 



with the function go defined in ([3|). For Na clusters, we 
have £f / Sum = 0.93 and our criterion yields a value Ae c ~ 
(3/5)e F . 

We have verified the robustness of our criterion ([2"6"j) 
by exploring different physical parameters, like the size of 
the cluster. Indeed, it can be seen in Figure [T] that above 
Ae c the RPA sum evaluated from (fT5|) (solid line) and 
the estimate (|2"4"|) (dashed line) are close to each other, 
showing that for high energies the RPA sum is essentially 
given by the contributions coming from the low-energy p- 
h excitations. We have checked that this feature of Ae c is 
independent of the size a of the nanoparticle. 



In order to have a well-defined collective excitation, the 
cutoff energy Ae c must obviously be larger than the min- 
imal p-h excitation energy (|14[) . For Na, we find that this 
condition is already verified with only N — 20 conduction 
electrons in the nanoparticle, in agreement with the ex- 
periments of reference [7] and the numerical calculations 
of reference [8] . 

The splitting in low and high-energy p-h excitations is 
important in order to justify the separation of the elec- 
tronic degrees of freedom into centre-of-mass and relative 
coordinates [5,6,13,18]. Within such an approach, the 
Hamiltonian of the electronic system is written as 



H = H r 



Hrel + He 



(31) 



The first term describes the centre-of-mass motion as a 
harmonic oscillator. The second term describes the rela- 
tive coordinates as independent fermions in the effective 
mean-field potential. The coupling term H c describes the 
creation or annihilation of a surface plasmon by destruc- 
tion or creation of p-h pairs. Its strength is given by the 
coupling constant A = (hm c ujf A /2N) 1 ^ 2 and the dipole 
matrix element d p h [6]. It is important to remark that, 
even if we have used the separability hypothesis of the 
residual interaction to justify the decomposition pip , the 
model that is put forward does not rely on such an ap- 
proximation. 



5 Dynamics of the relative-coordinate system 

The decomposition (j3"Tj) of H suggests to treat the collec- 
tive coordinate as a simple system of one degree of freedom 
which is coupled to an environment with many degrees of 
freedom. The latter are the relative coordinates described 
by Hrei. This is the approach taken in [20]. In this pic- 
ture, the time evolution of the centre-of-mass system (i.e., 
the surface plasmon) strongly depends on the dynamics of 
the relative-coordinate system. Such a dynamics is char- 
acterised by a correlation function which can be written 
at zero temperature as [20] 

C(t) =A 2 J2 \d ph \ 2 e iA ^ h O(Ae ph - Ae c ) . (32) 
ph 

In order to evaluate the above expression, it is helpful to 
introduce its Fourier transform 

E(Ae) = ^-A 2 \d P h\ 2 S(Ae - Ae ph )0(Ae ph - Ae c ) 

ph 

(33) 

which has been calculated in [6] . Finite temperatures were 
shown to result in a small quadratic correction. Consis- 
tently with the results of the preceding sections, we em- 
ploy our low energy estimates (f2Tj) and ([22]) to obtain 



(34) 



This result is consistent with the one of reference [6] (see 
Eq. (34) in there) in the limit Ae -C fr^M and for zero 
temperature. 
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In principle we could calculate C(t) by taking the in- 
verse Fourier transform of (|34|) . The decay of C(t) for very 
long times is dominated by the discontinuity of S at Ae c . 
This is somehow problematic since Ae c can only be es- 
timated as we did in Section 31 and since the functional 
form of the long time decay depends on how sharply the 
cutoff is implemented. However, it is important to realise 
that it is not the very long time behaviour that determines 
the relevant decay of the correlation function, but rather 
the typical values at which C(t) is reduced by an impor- 
tant factor from its initial value C(0). We then estimate 
the correlation time as the mean decay time of C (t) , 



d fC(t) 



d "d* vno, 
1 



C(0) 



dtC{t) 



i) 



(35) 



If C(t) were an exponentially decreasing function, (r cor ) 
would simply reduce to the inverse of the decay rate. Using 
the definition (|3"3"|) . we get 



/>oo 

h / dAeS{Ae)/Ae 

JAe c 

/>oo 

/ dAeS(Ae) 



(36) 



Given the fast decay of the function S ~ l/Ae 3 , the above 
integrals are dominated by their lower limit Ae c and we 
have 

<*->■»§£. <m 

Since Ae c is of the order of (3/5)ef, we see that the re- 
sponse time (or correlation time) of the electronic envi- 
ronment is of the order of its inverse Fermi energy. 

The estimation of the characteristic response time of 
the electronic environment is crucial in justifying the Mar- 
kovian approximation used in [20]. In that work, the de- 
grees of freedom corresponding to the relative coordinates 
were integrated out and treated as an incoherent heat bath 
that acts on the collective coordinate. Such an approach 
relies on the fast response of the environment as com- 
pared to the time evolution of the surface plasmon. The 
typical scale for the latter is the inverse of the decay rate, 
T sp = I/7. Using {2J and (f30|) . we therefore have 



' sp 



kid 



M 



9[.9o(ef/^m)] 2 £f 



k-pa . 



(38) 



For the example of Na nanoparticles worked in Section [3l 
we have r S p/(r cor ) ~ fcpa. This is a safe limit since in not 
too small nanoparticles, fcpa ^> 1. As the size of the cluster 
increases, the applicability of the Markovian approxima- 
tion is more justified. This is expected since the electronic 
bath has more and more degrees of freedom, approach- 
ing an "environment" in the sense of quantum dissipation. 
Since the physical parameters of alkaline nanoparticles en- 
tering (|38|) are close to that of noble-metal clusters, the 
dynamics of the surface plasmon can be expected to be 
Markovian in that case too. 



6 Conclusions 

We have studied the role of particle-hole excitations on 
the dynamics of the surface plasmon. A key concept in 
this analysis is the separation into low-energy excitations 
which lead to the collective excitation once they are mixed 
by the residual interaction, and high-energy excitations 
that act as an environment damping the resonance. Using 
the random phase approximation and assuming the sep- 
arability of the residual interaction, we have established 
a criterion for estimating the cutoff energy separating the 
low- and high-energy subspaces. The resulting cutoff en- 
ergy is approximately (3/5)ef for the case of Na nanopar- 
ticles. 

Since the number of electrons in the cluster is finite, 
the assumption that the high-energy particle-hole excita- 
tions act on the collective excitation as an environment, 
introducing friction in its dynamics, may be questionable. 
What settles this issue is the ratio between the typical evo- 
lution time of the collective excitation and the one of the 
high-energy particle-hole excitations. The former is given 
by the inverse of the plasmon linewidth, while the latter 
is obtained from the decay of the correlation function of 
the environment. We have found that this ratio improves 
with increasing cluster size. Even for a small cluster with 
a = 1 nm, the ratio is approximately 10, justifying the 
use of the Markovian approximation which assumes a fast 
time evolution of the environment with respect to the one 
of the collective excitation. 

The relevance of memory effects in the electronic dy- 
namics of small clusters is of current interest, due to the 
advance in time-resolved experimental techniques [14-17]. 
First-principle calculations have recently addressed this is- 
sue by comparing time-dependent density functional the- 
ories with and without memory effects for small gold clus- 
ters [25]. For very small clusters (N < 8) memory effects 
were shown to be important. We stress that the mem- 
ory considered in [25] is that of the electron gas as a 
whole, while we are concerned in this work with the mem- 
ory arising from the dynamical evolution of the relative- 
coordinate subsystem. It would be interesting to consider 
cluster sizes intermediate between the ones considered in 
the present work and those of reference [25] in order to 
study the emergence of memory effects. 
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A Lowest energy of the particle-hole 
spectrum 

If we consider the cluster as a hard- wall sphere of radius a, 
its eigenstates are given in terms of spherical Bessel func- 
tions. Using the large ka expansion of the latter (semiclas- 
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sical high-energy limit), the quantisation condition reads 



ka = 7T 



(39) 



with I and n non-negative integers. The energy of a single- 
particle (hole) state is related to its wavevector k p ^) , its 
total angular momentum l p (h), and its radial quantum 



number n 



■ P (h) 



as 



£ p(h) ~ 



h2k p(h) 

2m c 



kpa 



l p(h) 



l P(h) 



(40) 



Thus, the energy of a p-h excitation entering the RPA sum 
JH is 



Ae p h — ef 



k-pa 



lp lh 

2 

l P + hi 



n p - n h 
+ n p + n h 



(41) 



Notice that using the exact quantum mechanical spec- 
trum in (fl3|) would not change significantly the result de- 
picted in Figure [TJ since the approximation (|39|) is very 
reliable for states close to the Fermi energy. We have also 
checked that generating the p-h excitation energies ran- 
domly in the RPA sum Q does not affect the physical 
picture of Figure Q] Indeed, the main ingredient to under- 
stand such a picture is the fast decay of the dipole matrix 
element with the p-h energy. 

The expressions for the dipole matrix element (TT2"|) as 
well as for the typical dipole matrix element (|2"2"|) diverge 
in the limit of a small p-h energy. It is therefore crucial 
for our analysis to determine the appropriate minimal p-h 
energy Z\e m m that renders this divergence unphysical. 

This can be achieved by imposing the dipole selection 
rules in (1411) and that the energy difference is minimal. The 
first condition dictates that lh = l p ± 1 and rrih — m p . 
Therefore there are two ways of obtaining the minimal 
energy difference: n p = with l p = l^ + l and n p = n^ + l 



with l p = lh — 1. In both cases we have 



kpa /it kp 



(42) 



Since we are interested in states close to the Fermi level, 
we can simplify (|42|) to expression (|14jl . 

If we consider sodium clusters with kpa = 30 (a = 
3.3 nm and N ~ 4000 conduction electrons per spin direc- 
tion), we have Z\e m i n ~ £f/10. This is a much larger en- 
ergy than the lowest one we can observe in the numerically 
generated excitation spectrum (see Fig. 1 in Ref. [12]). 
However, the two results are reconciled once we take into 
account the large degeneracy yielded by our approximate 
quantisation condition (I39p . 

The degeneracy of p-h excitations with minimal energy 
is given by twice the number of pairs (lh,rih) compatible 



with kh = k-p and Ae p h = Ae m i n . Indeed, we have seen 
that there are two possible particle states p starting from h 
and verifying the above-mentioned conditions. For each n 
between 1 and kpa/n, there is a value of I = 2{kpa/-K — n) 
and therefore the number of degenerate p-h excitations 
with energy Z\£ m in is 

kpa/n 

JV = 2 (2/ + 1) 

n=0 
kpa / 7r 

= 2£ 



7T 



kpa 



(43) 



This degeneracy factor has to be included in Figure [TJ 
and it is crucial for the determination of the collective 
excitation. 



B Local density of the dipole matrix element 

Equation (fTH|) defines the local density of dipole matrix 
elements connecting states at energies e and e + As. We 
used particle and hole states in our definition, since this 
is the main interest of our work. But note that the calcu- 
lation presented in this appendix is not restricted to that 
case and can be easily extended to any states. The result 
would be of course unchanged. 

Local densities of matrix elements of arbitrary oper- 
ators have been thoroughly studied as they can be eas- 
ily connected with physical properties, ranging from far- 
infrared absorption in small particles [26] to electronic 
lifetimes of quantum dots [27] . A semiclassical theory for 
the local density of matrix elements has been developed 
[28-30], where (|16p can be expressed as a smooth part 
given by correlations along classical trajectories plus a pe- 
riodic orbit expansion. We will not follow here this general 
procedure, but use the simple form of the dipole matrix 
elements (fT0|) for states confined in a hard-wall sphere 
and the semiclassical approximation applied to the radial 
(fixed I) problem [12,13]. 

Introducing the Z-fixed density of states, which in lead- 
ing order in h is given by 



Qi{s) 



^2m c a 2 e/h 2 - (I + 1/2) 2 
2tt£ 



(44) 



we can write with the help of equations (fTT)l) - (fT2"|) 

C(e,Ae) 



2h 2 \ 2 ee' 



m c a J 3Ae 4 

x Qh( £ ) [{h + l)£^+i(e') + lhQl h -i{ £ ')] 
i h =o 

(45) 

where Z max is the maximum allowed lh for an energy e, 
while e 1 — e + As. In the semiclassical limit we can take 
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//,, ~ lh + 1 3> 1 and convert the sum into an integral. 
Thus, 



, , 1 [2h 2 Y 1 



5— e /2 

an 



2ma 2 



(46) 



Performing the remaining integral over the angular mo- 
mentum I finally yields the result (fTTjl . 



C Density of particle-hole excitations 

The density of p-h excitations with energy Ae is defined 
in (T2I7|) and can be written as 

P p - h (Ae) = f ¥ de h Y,{Uh + l)& h (e h ) 

Je^ — Ae 1 

x [gi h+1 (e h + Ae) + gi h -i(e h + As)] . (47) 

Using the semiclassical density of states (j4"4")) and perform- 
ing the sum in the limit 1^ ^> 1, we obtain 



p v- l \Ae) 



Ae 2 [2m c a 2 \ 2 f £F 



8ir 2 V h 2 



de h 



F(e h /Ae) 



-Ae e h {e h +Ae) 



(48) 

where the function F has been defined in (|18p . Performing 
the remaining integral over the hole energy in the limit 
Ae <C £f is straightforward and leads to the result (|2Tj) . 
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